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I. Introduction 



Numerical simulations of the freely decaying, incompressible Navier-Stokes equations in two 
dimensions have shown that under appropriate conditions and after a relatively short period of 
chaotic mixing, the vorticity becomes strongly, localized in a collection of vortices which move 
in a background of weak vorticity gradients p\. As long as their sizes are much smaller than 
the extension of the domain, the collection of vortices may evolve self-similarly in time P ,i ,@ 
Ji J§ j3] until one large-scale structure remains. If the corresponding Reynolds number is large 
enough, the time evolution of these so-called coherent structures is usually given by a uniform 
translation or rotation and by relatively slow decay and diffusion, the last two are due to the 
presence of a non-vanishing viscosity. In other words, in a co-translating or co-rotating frame 
of reference, one has quasi-stationary structures (QSS) which are, to a good approximation, 
stationary solutions of the inviscid Euler equations. Accordingly, their corresponding vortic- 
ity^ fields ujs{x,y) and stream functions ips{,x,y) are, to a good approximation, functionally 
related, i.e., u!s{x,y) ~ ujs{ips{x,y)). Similar phenomena have been observed in the quasi two- 
dimensional flows studied in the laboratory P The only exception to this rule is provided 
by the large-scale, oscillatory states that occasionally result at the end of the chaotic mixing 
period many cases, e.g., when the initial vorticity field is randomly distributed in 

space, the formation of the QSS corresponds to the segregation of different- sign vorticity and 
the subsequent coalescence of equal-sign vorticity, i.e., to a spatial demixing of vorticity. 
Besides the theoretical fluid-dynamics context, a good understanding of the above- described 
process has implications in many other physically interesting situations like: geophysical flows 
f^, plasmas in magnetic fields ^ ], galaxy structure [0], etc. For these reasons numerical 
and experimental studies are still being performed and have already led to a number of "scat- 
ter plots", i.e., to the determination of the ws-ips functional relation as a characterization of 
the QSS which appear under different circumstances. Simultaneously, on the theoretical side, 
approaches have been proposed which attempt at, among other things, predicting the QSS 
directly from the initial vorticity field; if successful in this, such methods would also alleviate 
the need of performing costly numerical and laboratory studies. 

The above-mentioned studies point out the large enstrophy decay that often takes place 
during the formation of the QSS; sometimes, also the evolution of the skewness is reported. 
But for these two lowest-order moments, little attention has been paid to the evolution of the 
vorticity-area distribution, defined in equation (|1|), during the formation of the QSS. In the 
context of 2D flows, this distribution plays a very important role: with appropriate boundary 
conditions, it is conserved by the inviscid Euler equations and the stationary solutions are the 
maximizers of the energy for the given vorticity distribution [111 ,0], see also Subsection [1II|A. 



In the present work we study the time evolution of the vorticity-area distrbution in two- 
dimensional, incompressible and viscous fluids. Many of the ideas we present should be appli- 
cable also to more realistic systems, e.g., when potential vorticity is the Lagrangian invariant. 
The paper is structured as follows: in the following Section, we derive, from the Navier-Stokes 
equation, the time evolution of the vorticity-area density. It is an advection-diffusion equation 
with a time dependent, negative diffusion coefficient. For the purpose of illustration, explicit 
calculations are presented in Subsection |11|B| for the case of a Gaussian monopole and, in Subsec- 
tion PCL for the case of self-similar decay described by Bartello and Warn in fl. Based on these 



^We use the subscript S in order to indicate that the field corresponds to an observed QSS. 
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diffusion coefficients, it would be interesting to set up a classification distinguishing different 
various typical scenarios leading to the QSS. Considering the QSS, a very natural but difficult 
question arises about its relation to theoretical predictions, namely, it is not trivial to quantify 
howgood or bad the agreement between observation and prediction is, as already stressed, e.g., 
in |0 ] . For this and related reasons, in Section |T|, we show that a perfect agreement between 
an observed QSS and the corresponding prediction obtained through the statistical-mechanics 
approach as developed by J. Miller et al. [0 ,0 ,0] and by Robert and Sommeria ^ ,0] would 
imply the equality of the difference in the moments of the initial and final vorticity distributions 
on the one hand and a set of quantities that can be directly obtained from the experimental 
uJs-ips relation on the other side. The details of the proof can be found in the Appendix. In 
II HQ we discuss how to use these quantities as yardsticks in order to quantify the validity of the 
statistical-mechanics approach in numerical and laboratory experiments. In the last Section we 
summarize our results and add some comments. 



II. Vorticity-area distribution 

A. Time evolution 

It turns out that the vorticity-area density undergoes an anti-diffusion process as we next 
show. The time dependent vorticity-area density G{a, t) is given by 

G{a,t) := / dxdy 5{a — Lu{x,y,t)), (1) 

J A 

where A denotes the domain and the vorticity field u{x,y,t):=k ■ (V x tJ) evolves according to 
the Navier- Stokes equation 

du ^ „ ^ , , 

— + vVuj = z/Acj, (2) 
ot 

where the incompressibility condition V-t? = has been taken into account. The Navier-Stokes 
equation determines the time evolution of the vorticity-area density; one has 

— — — = / dxdy 5'{a — iu{x,y,t))v ■ Viu — I' / dxdy 5'{a — iu{x,y,t))Auj 

J A J A 

f ^ d f 

= — dxdy V ■ VSi^cr — Lj{x,y,t)) — u— dxdy 6{(r — Lj{x,y,t))Au. 

J A C"^ J A 

We will assume impermeable boundaries dA, i.e., that the velocity component perpendicular 
to dA vanishes. Therefore, the first integral in the last expression is zero. Partial integration 
of the second integral leads to 

dG(a,t) d"^ [ 1 1 ^/ , a 9 f ^, ^ „ 

— — z/7— — dxdy d[(j — uj[x,y,t)) \Vuj\ — i^t^ f o{a — ujvuj ■ ndl. [6) 



dt da'^ Ja ' ' da 



The last term represents the net vorticity generation or destruction that occurs on the boundary 
OA of the domain, n is the outward oriented, unit vector normal to the boundary. This source 
term vanishes in some special cases like doubly periodic boundary conditions or when the 
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support of the vorticity field remains always away from the boundary. 

From the definition of G{a,t) one sees that if G{a,t) = and \Vuj\ is finite, then also the 
integrals in the last expression must vanish. Consequently, we can define an effective diffusion 
coefficient D{a,t) through 

z/ dxdy6{a -uj{x,y,t))\Vuj\'^ =: D{a,t)G{a,t), (4) 

J A 

From this definiton it is clear that D{a, t) > is the average of z/ | Va;|^ over the area on which 
the vorticity takes the value a. Therefore, the time evolution of G{a, t) can be written as an 
advection-diffusion equation in cr-space with a source term, i.e.. 



dG{a,t) _ d 
di ~ ~da 



s{a, t) + —^^G{a, t) + D{a, t) ' 



(5) 



with a negative diffusion coefficient —D{a,t), a "velocity field" dD{a,t)/da and where minus 
the cr-derivative of s(cr, t) := u jg^ 6{a — u!)Vuj ■ ndl is the vorticity source at the boundary. 
Introducing the vorticity moments 



r„(t) := jdaa'^Gia,t) = j dxdy uj'^{x,y,t), 



one can check that the equations above imply that the first moment of the distribution G{a, t), 
i.e., the total circulation Fi, evolves according to 



dVi 



da s{a, t), 



dt 

and that the even moments T2n{t) '■= J da a'^"'G{a,t) change in time according to 



dT 



dt 



2n j da a'^-^sia, t) - u2n{2n - I) j dxdy u^^''-^^ | VcjI^ 



In particular, when the boundary source term s{a,t) vanishes, the total circulation Fi is con- 
served and the even moments decay in time 



dT 



2n 



dt 



-v2n{2n - I) j dxdy u'^^'^'^^ \Vuj\^ < 0. (6) 



Therefore, in cr-space one has anti-diffusion at all times and with \Wuj\ — *>t^oo when the 
boundary source term s{a,t) vanishes, the final vorticity-area distribution is 

lim G{a^ t) = A6{a — Co) and lim D{a, t) = 0, 

t~*oo t^oo 

where A is the area of the domain and ui is the average vorticity, i.e., uj = Ti/A. 

In most cases, it is not possible to make an a priori calculation of the effective diffusion 
coefficient D{a,t). On the other hand, its computation is straightforward if a solution of the 
Navier-Stokes equation is known. These diffusion coefficients and, in particular the product 
D{a,t)G{a,t), confer Figs. 1 and 2 in |y, may lead to a classification of different scenarios for 
and stages in the formation of QSS. It is also worthwhile recalling that conditional averages very 
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similar to D{(T,t), confer (H), play an important role in, e.g., the advection of passive scalars 
by a random velocity field. In the case of passive-scalar advection by self-similar, stationary 
turbulent flows, Kraichnan has proposed a way of computing such quantities P]. 
For the purpose of illustration, in the next Subsection, we use an exact analytic solution in 
order to compute the corresponding G{(J, t) and D{a, t). Another tractable case occurs when the 
flow evolves self-similarly; in Subsection |lip| we apply these ideas to the self-similar evolution 
data obtained by Bartello and Warn [0]. 

B. A simple example 

As a simple example consider the exact solution of the Navier-Stokes equation in an infinite 
domain given by a Gaussian monopole with circulation Fi, 

u!g{x, y, t) = Fi {Anut)^'^ exp(— r^/4i^t), with := x'^ + and t > 0. 

Then, in cylindrical coordinates (r, (/)), 



S{a — ujg{x, y, t)) dxdy = 5(r^ — -R^(cr, t)) -dr'^dcj)/ 



dujG 



1 

2/ -x . .. ,(yAT^vt, 



where i^G{x,y,t\^j^ = a, i.e., R {a,t) := -4z/tln(— - — ) 

J- 1 



duG 

and 



a 



and we have that for such a Gaussian monopole, the vorticity-area density is 

GG{(y,t) = 0, for a < 0, 

= 4:7iuta~^, for < cr < crmaxit), 
= 0, for amax(i) < cr, 
with CTmax(t) = ^1 (47rz/t)"^ . 

The divergence at o" — ^ is due to the increasingly large areas occupied by vanishingly small 
vorticity associated with the tails of the Gaussian profile. Due to this divergence, the density 
GG^cr, t) is not integrable as it should since the domain A is infinite. In spite of this divergence, 
all the (T-moments are finite, in particular, the first moment, i.e., the circulation, equals Fi and 
the second moment, i.e., the enstrophy, is F^/Svrz/t. As expected, the circulation is constant in 
time while the enstrophy decays to zero[|. It is also interesting to notice that while the maximum 
vorticity value, crmax(^) = Ti/inut, occupies only one point, i.e., a set of zero dimension, 
the density GG{cr,t) remains finite for a o"max(^), more precisely, limo.^o.^^^(() Gg(o", t) = 
(47ri/t)Vri. 

Moreover, in this simple example one can also compute 



V I dxdyS{a — ujg{x, y, t)) IVc^gI^ = 47rz/cr In ( ^ ^ ) , for < a < crmax(^), 



^max (^) 



■^By contrast, the spatial second moments of Lu{x,y,t) increase in time like 2i/t, e.g., / ^ dxdyx^ujG{x,y,t) 
2vt. 
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Figure 1: Plot of the dimensionless quantities (upper curve) {4:nh'amax)~^Dc{(T,t)GG{c,t) and 
(lower curve) Ti^ini'o-f^^^)'^ DG{cr,t) , both as function of x = o"/(7max in the case of a Gaussian 
monopole. 



so that the corresponding effective diffusion coefficient is 

DgM = y In (^^I^^ , for < a < a^^t)- 

The vanishing of this Dg{(J, t) with a corresponds to the vanishingly small spatial gradients 
of vorticity at large distances from the vortex core; this gradient vanishes also at the center of 
the vortex leading to a (weaker) vanishing of Dg{(t, t) for a y CTma^it) = Ti/Anut. The Gaussian 
vortex is totally dominated by viscosity and yet the corresponding effective diffusion coefficient 
DgIct, t) is not proportional to the viscosity v, as one would naively expect, but to In u. In Fig. 1 
we plot the dimensionless quantities (47rz/crmax)~^ -Dg(o", t)GG{o', t) and Fi {Airua'^^^) ^ -DG(cr, t) 
as functions of o ja^^^. 

Fig.l. Plot of the dimensionless quantities (full curve) (47rz/cTmax)^^-DG(c", t)GG(o"; ^) and 
(dotted curve) Fi(47rz/(T^g^^)~^Z)G(C) ^)) both as function of a; = a/a^^ in the case of a Gaussian 
monopole. 



C. Self-similar decay 

In the case of a collection of vortices that evolves self-similarly in time, it is convenient to 
introduce the dimensionless independent variable ^ := {a — uj)t and the dimensionless functions^ 
G(^) := G{a — LD,t)/At and D{^) := D{a — uj,t)/ {a — lD)'^. When the boundary source term 
s((j, t) vanishes, as it does in the case of a doubly periodic domain or when the vorticity support 
remains well separated from the boundary, the self-similar form of (^ is 

or 

^We assume that A is finite and that t > 0. 



7 



Assuming that there are no singularities in the vorticity field, one has t) ^ and it 

|(t|— >oo 

follows then that the constant in the last expression must be zero. Measuring the self-similar 
density G{C,), one can solve the last equation for the corresponding diffusion coefficient D{^) 
and get 



1 



D(n = . / dssG(s), 

where the value of the lower limit of integration b must be chosen according to an appropriate 
"boundary condition" as illustrated below. If for large |^| the dimensionless vorticity density 
G(^) decays algebraically like |^|~^" (see below) and 2a < 3, then it follows that the most 
general decay of D{^) is of the form |^|~^ + a |^|^^"~^'* with a an appropriate constant. 

We apply these results to two specific cases: 1) If the self-similar vorticity distribution 
happens to be Gaussian, i.e., if G{^) = (^oV^27r) exp (— ^^/2^^) . Then, with D(±oo) = 
as "boundary condition" one obtains that -D(0 = ^o^~^ or, going back to the original quan- 
tities, G{a,t) = At{^oV27i)' exp [— t^((T — Ci;)^/2^o] and the negative diffusion coefficient is 
—D{cx,t) = — ^ot^^. This is the only self-similar case with a cr-independent D {a,t) . In this 
case, the time-evolution equation takes a particularly simple form, namely 

dG{a,t) ^ d'G{a,t) 

dt ^° da^ ' 

. dG{a, r) 1 d^G(a, r) . , , 

In agreement with our findings in IIA, we see that in this case the squared width of the Gaussian 
distribution decreases in time like r = 

2) The second application is to the self-similar distributions found by Bartello and Warn 
in their simulations performed in a doubly periodic domain of size A. Qualitatively speaking, 
their results can be summarized by the following expression 

Gsio = c{eo+er\ \^\<^m (7) 

GsiO = 0, > ^M. 

with ~ 10, a ^ 0.7 and growing^ approximately like ^/i, from 200 to 500. The value 
of c being such that J dC, G{^) = 1. Vorticity values such that |cr| t < C,o are associated mainly 
with thin filaments in the background "sea" while those such that |cr| t > C,o correspond to the 
localized vortices. At the positions with the largest vorticity value the gradient Vuj vanishes, 
therefore, it is natural to take as "boundary condition" Ds{C,m) = 0. One gets then the following 
effective diffusion coefficient in vorticity-space 

Going back to the original variables, this effective diffusion coefficient reads. 



^This time-dependence destroys the exact self-similarity of these solutions. 
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Figure 2: Plot of the dimensionless quantities (full curve) t^Ds(cr, t)^^/ and (dotted curve) 
4 ^cA^j^^~"''^o°'j t^Ds{a, t)Gs{a, t) as functions of x = at/^M in the case of self-similar decay 
with S,o = 10 and for S,m = 300, as discussed in the main text. 



The average of IVc^l in the thin filaments is not zero so that at cr = we have 

So 



D,(0,t) 



2(l-a)t3 



1 + 



C2 \ 
SM 



So 



1 



f2af2(l-a) 
So SM 



2(l-a)t3- 



The origin, o" = is a local minimum of Ds(cr, t), moreover there is one maximum, namely 



max D(a,t) = a^-^- 



2t3 



at a*t=JaT=^ieo+eMit))-eo^c^^i 



For large t, the maximum decays like t ^ while D{0,t) decays faster, namely like t 



M- 



(2+a) 



In Fig. 2, we plot the dimensionless quantities t^Ds{a, t)^j^,/ and 4 ( cA^'^^j^ "^C.o" ) t^Ds{a, t)Gs{cr, t) 
as functions of x = at/^M- 



III. The changes in the vorticity-area distribution 

A. General considerations 

In comparison to the case of a passive scalar, the spatial distribution of vorticity does not 
play such a central role as one realizes from Arnold's observation that the stationary solutions 
of the Euler equations in two dimensions correspond to energy extrema under the constraint 
of fixed vorticity areas [El ,0]. Arnold's observation says then that the stationary solutions of 
the 2D Euler equations, ujs{x, y), are the states with extremal values of the energy compatible 
with the vorticity density 




dxdy5{a — ujs{x, y)). 
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Therefore, the vorticity-area density Gs{(t) of a QSS (and the geometry of the domain) de- 
termine uJs{x,y), the spatial distribution of vorticity in the coherent structure. From this we 
conclude that when studying the process leading to the QSS in viscous fluids, special atten- 
tion should be paid to the differences between the initial vorticity-area density G{a, 0) and 
Gs{cf), the one in the QSS. A convenient way of studying these changes would be through the 
differences in the moments of these distributions, which will be denoted by A„, i.e., 

A„, := - rf , with 



r° := 

n 



r 



dxdyuj"^{x,y,0) ^ J da a'^G {a, 0) and 
dxdyuj'^{x,y) = da a'' G sic). 



The dimensionless ratios A„/r^ would offer a good characterization of the changes experienced 
by the vorticity-area distribution. In the next Subsections we present another possibility, linked 
to the predictions of a QSS according to a statistical-mechanics approach, which is constructed 
from an cusH^s) relation measured either in experiments or in numerical simulations. 

B. The changes in the moments according to the statistical mechan- 
ics approach 

It is proven in the Appendix that when the quasi-stationary vorticity field ujs{x,y) and 
the initial field u{x, y, 0) are related as predicted by the statistical mechanical theory then the 
observed A„ take the values 5„ 



5„ = y" dxdyiniip), 
with ii := and 



(8) 



1^' 



defined in terms of the associted VL{%1)) relation and an inverse temperature j3. In particular, for 
n < 4, we have 



5i = 0, 

§2 = -13-^ jdxdy {dVt/dip), 



(9) 



S3 — P ^ J dxdy 
84 — —(3'^ I dxdy 



d^n 3/50^^' 
dip'^ dip 



In the next Subsection we propose the use of 5n as yardsticks in order to quantify the departure 
of the observed changes A„ from the theoretical predictions. 
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C. Analysis of experimental results 

In many cases, the predictions obtained from different theories are not very different and 
it is not obvious which prediction agrees better with the experimental data, see, e.g., [0 JhI] . 
Therefore, it is important to develop objective, quantitative measures of such an agreement. 

The results of the previous Subsection lead us to conclude that there is useful information 
encoded in the functional dependence of uJsi^^s) and that this information can be used for the 
quantification of the vorticity redistribution process in any experiment or numerical simulation, 
i.e., also when (|1^) and do not necessarily hold, as long as there is no leakage and creation 



or destruction of vorticity at the boundary is negligeable. We propose that one should proceed 
as follows: 

1) Identify the predicted fi(V') relation of the preceeding Subsection and the Appendix with 
the ujs{x,y) ~ ujsi'ips) of the observed QSS; usually, this satisfies dus/dips 7^ 0, 

2) Determine an effective value of (3 from A2, the measured change in the second moment of 
the vorticity-area distribution and from the measured uJs{ips) relation by, confer the second 
line of (I), 

^ jdxdy{dujsldii)s) ^^^^ 
Since A2 > 0, the sign of (3 is always opposite to that of dusjdips- 

3) Using this value of (3 compute from equation (H), with VLi^ip) replaced by ujs{jps)i the values 
of the yardsticks 5^ for n > 3. 

4) The measured changes in the third and higher moments, A„ with n > 3, should be quantified 
by the dimensionless numbers an '■= A„/(5„. These numbers are all equal to 1 if and only if the 
QSS agrees with the statistical mechanical prediction corresponding to the initial distribution 
G(a, 0) with equal initial and final energies, 

^s^i^s) corresponds to a statitistical equilibrium <-^a3 = -- - = a„ = -- - = l. (11) 

It is for this reason that one should prefer these dimensionless quantities to other ones like, 
e.g., A„/rf. 

In closing, it may be worthwhile recalling that, at least in some cases, the agreement between 
an experiment and the statistical mechanical prediction can be greatly improved by taking 
as "initial condition" not the field at the start of the experiment but a later one, after some 
preliminary mixing has taken place but well before the QSS appears, see pi]. This improvement 
can be quantified by measuring the convergence of the corresponding a„ = A„/(5„ towards 1. In 
other detailed consideration of the boundary is necessary and, sometimes, the satistical 

mechanics approach may be applicable in a well-chosen subdomain, see 

D. Examples 

A possible Vl{'ip) relation resulting from the statistical mechanical theory that may be com- 
pared successfully to many experimentally found curves is 



(12) 
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with appropriately chosen constants f^o, B, C and x- In particular, the flattening in the scatter 
plots observed in pzt] can be fitted by this expression while the case C = 0, corresponds to 
the identical point- vortices model pi ,0] and the case x ^ ^ with x^o/ {B + C) —infinite, 
corresponds to a linear scatter-plot. In all these cases, the 5„ can be derived, as explained in 
the Appendix, by means of a cumulant generating function, confer ([T7|), 



pQo, B + Ccoshx{ip + X) 
— In ■ 



xC B + C cosh xi' 

Expanding this in powers of A leads to the cumulants of the microscopic vorticity distribution. 
In particular, 

xflo B cosh xi^ + C 
13 (5 + Ccoshx^)^' 

Integrating this over the area, one obtains that 

X^o f J J B cosh x'ip + C 

\2- 



^2 = j ^^^y 



{B + C coshx^py 

Knowledge of and of fi* ('?/') allows us to determine the value of /?. Once /? has been fixed. 



one can apply (D) and ([IT]) in order to quantify the higher-order moments and in so doing to 
estimate the validity of equation (|12|) as the prediction from the statistical-mechanics approach. 
If the scatter-plot is linear, = kip, then (O) tells us that 



This implies a simple relation between the initial vorticity-area distribution G{a, 0) and Gs{(j), 
the distribution in the QSS, namely 

J da exp{-X/3a)G (a, 0) = exp{^X'^pk) j da exp{-Xpa)Gs (a) . 

Recall that, as already indicated immediately after (p!OD, pk < 0. Assuming that the Laplace 
transformation can be inverted, we get for the linear case = kip, that 

G(a,0)= JdT exp(-i^^)G5(r) 

In particular, if the QSS has a Gaussian vorticity-area density with a variance equal to and 
the predictions of the statistical mechanical approach hold, then the initial distribution G {a, 0) 
must also be a Gaussian with variance equal to (S^ - (3k) > T?, in agreement with the results 
m Notice that knowing G{a, 0) and the energy of the initial state is not enough in order 

to determine the spatial dependence of the initial vorticity field u!{x,y, 0). 
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IV. Conclusions 

In this paper, the object at the center of our attention has been the vorticity-area density 
G{a,t) and its time evolution in two-dimensional, viscous flows. In Section |I| we have shown 
that this density evolves according to an advection-diffusion equation, equation (^, with a 
time dependent, negative diffusion coefficient. If vorticity is destroyed or created at the domain 
boundaries then the evolution equation contains also a source term. The equation is exact: 
it follows from the Navier-Stokes equation with no approximations made. For the purpose of 
illustration, explicit calculations have been presented for the case of a Gaussian monopole in 
and for the case of self-similar decay in |11|C| . We think that it will be instructive to apply 



these ideas to the analysis of data that is available from numerical simulations and laboratory 
experiments. In fact, it should be possible to determine this effective diffusion coefficient -D(cr, t) 
on the basis of such data. Then it would be of interest to establish a quantitative classification 
of the QSS formation processes, e.g., by considering the various possible behaviours of this 
coefficient and, in particular, confer Figs. 1 and 2, of G{(j,t)D{a,t). In the case of self-similar 
decay, one could attempt a closure approximation in order to predict the effective diffusion 
coefficient like it is done, for very similar quantities, in the theory of passive-scalar dispersion 
by random velocity fields, e.g., Kraichnan's linear Ansatz pi]. 

In Section |T| we considered the changes in G{a, t) when starting from an arbitrary vorticity 
field and ending at a high- Reynolds' number, quasi-stationary state characterized by an usi^ps) 
relation. In |111|B| and in the Appendix, we showed how to generate from such an u-ip plot 



an infinite set of moments 5„, confer (P). The changes A„ in the moments of the vorticity 
distribution that are observed in a numerical simulation or in the laboratory equal these 6n 
if and only if the initial and final distributions are related to each other in the way predicted 
by the statistical mechanical approach. Therefore, these changes in the vorticity distribution 
moments can be quantified in terms of the dimensionless ratios A„/5„. The deviations of the 
ratios A„/5„ from the value 1 as determined on the basis of the data gives a direct way of 
quantifying the validity of the statistical-mechanics approach. In |III|C| we discussed how to 



apply this to experimental measurements provided that the leakage or creation and destruction 
of vorticity at the boundaries is negligeable. Finally, in Subsection |llip| , we considered two 



relevant u-ip relations: a linear one and a much more general one, given by equation (0). Many 
of these ideas should be applicable also to more realistic systems, e.g., when potential vorticity 
is a Lagrangian invariant. 
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from P.-H. Chavanis' comments and from the comments of the referees. The interest shown by 
H. Clercx, G.-J. van Heijst, S. Maassen and J.J. Rasmussen has been most stimulating. 



V. Appendix 

We prove now that if and only if the QSS happens to coincide with the prediction of the 
statistical mechanical approach, then the changes in the moments A„ take the values Sn as 
defined in (|]). 

Recall that in the statistical mechanical approach one identifies uJs{x,y), the vorticity of the 
QSS, with (cr) := fdaap{a,ilj), the average value of the microscopic vorticity a with respect 
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to a vorticity distribution p((j, ip) which is given by 

p{a, ip) := Z'^ exp [-[3ail){x, y) + , (13) 

with Z{ip) := J da exp [—Paip + 

and define ^{ip) := (a) = J daap{a,'ip). 

In (pTSD, /3 and /i(cr) are Lagrange multiphers such that the energy per unit mass and the 
microscopic-vorticity area distribution g{a) := j dxdy p{a,ip{x,y)) have the same values as in 
the initial distribution, i.e., g{a) = G{a,0). Consequently, the spatially integrated moments 
are given by 

j dxdy (a") = J daa''g{a) = J daa"Gia,0) = 1° (14) 

while ^n = j dxdy or {i)) = J dxdy (a)". 
Denote by 5„ the predicted change in the n-th moment, i.e., 

6n = I dxdy [{a"^) - {oY'] =: j dxdy in- (15) 



To the probability distribution p{(t,iIj), defined in (13), we associate a cumulant generating 
function 

k(A, ^) := In (exp (-A/3a)) . (16) 

This satisfies 

^^^ = -W + A). (17) 

as it can be shown by first noticing that 

(exp(-A/3a)) = ^^ 
so that fi;(A, '0) = lnZ{ip + A) — \nZ{ip) and then using J da ap{a,ip + A) = VL{tIj + A) for the 



computation of dn/dX. Expanding both sides of (|17D in powers of A, it follows that K,n{ip), the 
n-th local cumulant of — /5cr, is related to fl^ip) by 

'^nW = -P^.nii;). (18) 
For example, for 1 < n < 4, these equalities read 

P{a)=pn, 



2 ■ 



dip' 
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For our purposes, it is essential to eliminate from these equations products like fl (cr^) and (cr^) 
because their integrals over the whole domain cannot be related to known quantities. In fact, 
the differences (5„ can be directly expressed in terms of and its derivatives, i.e., in terms 
of known quantities, as we now show. To this end, one considers first the generating function 
of the local moment differences in, which will be denoted by i{X,ilj). One has that, confer (p!5|), 

z(A, ip) := (exp (— A/5(j)) — exp (—A/5 {a)) , 

This is related to the cumulants generating function /t(A,^/'), confer equation (0), by 

«(A, ij) = In [^(A, ij) + exp (-A/5 (a))] . 



Expanding both sides of this identity in powers of A and making use of (0), one obtains the 
recursive expressions for the in{ip) given in equation integrating these over the area A, one 
finally gets the results stated in [1I1|B . 
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